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Abstract: Measurements of the arrival directions of cosmic rays iiave not revealed their sources. High 
energy neutrino telescopes attempt to resolve the problem by detecting neutrinos whose directions are not 
scrambled by magnetic fields. The key issue is whether the neutrino flux produced in cosmic ray accel- 
erators is detectable. It is believed that the answer is affirmative, both for the galactic and extragalactic 
sources, provided the detector has kilometer-scale dimensions. We revisit the case for kilometer-scale 
neutrino detectors in a model-independent way by focussing on the energetics of the sources. The real 
breakthrough though has not been on the theory but on the technology front: the considerable technical 
hurdles to build such detectors have been overcome. 

Where extragalactic cosmic rays are concerned an alternative method to probe the accelerators consists in 
studying the arrival directions of neutrinos produced in interactions with the microwave background near 
the source, i.e. within a GZK radius. Their flux is calculable within large ambiguities but, in any case, 
low. It is therefore likely that detectors that are larger yet by several orders of magnitudes are required. 
These exploit novel techniques, such as detecting the secondary radiation at radio wavelengths emitted 
by neutrino induced showers. 



Cosmic Rays and Neutrinos 

An illustration of the neutrino sky is shown in 
Fig. 1 displaying a spectrum ranging from mi- 
crowave neutrinos produced in the big bang to 
GZK neutrinos associated with the highest energy 
cosmic rays[l]. The GZK neutrinos are the de- 
cay products of pions produced in the interaction 
of cosmic rays with microwave photons. These 
are the same interactions that shape the Greissen- 
Zatsepin-Kuzmin absorption feature in the spec- 
trum, hence their name. Prominently displayed 
in the figure is the flux of the highest energy at- 
mospheric neutrinos observed up to ~ 100 TeV 
by the AMANDA experiment [2]. This beam, 
very successfully mined for particle physics by the 
Superkamiokande-generation of experiments, will 
be exploited at yet higher energies [3]. Because of 
its steep spectrum, events above several hundreds 
of TeV become very rare, leaving a clear neutrino 
sky to be explored for sources of cosmic neutri- 
nos beyond the sun. Neutrino telescopes will open 
some ten orders of magnitude in neutrino wave- 
length, from their tens of GeV threshold to the EeV 



energy of GZK neutrinos.' The existence of cos- 
mic neutrinos with yet higher energy are a matter 
of speculation; they could be the decay products of 
cosmic remnants or topological defects associated 
with phase transitions in the early universe[4]. 

The neutrino fluxes anticipated [4] from non- 
thermal astronomical sources and, therefore, can- 
didate cosmic ray accelerators such as active galax- 
ies (AGN) and gamma ray bursts (GRB), dominate 
the atmospheric flux above ~ IPeV; see Fig. 1. 
The energy fluxes are at the level of the flux, as far 
as we know imobservable, of big bang neutrinos. 
The venture can nevertheless succeed by exploiting 
the relatively large high energy neutrino interaction 
cross sections in combination with detectors of gi- 
gaton size. Standard model physics is sufficient to 
establish that the cosmic fluxes shown are observ- 
able in a volume of 1 kilometer cubed instrumented 
with photomultipUers[4]. Neutrino telescopes de- 
tect the Cherenkov radiation from secondary par- 
ticles produced in the interactions of high energy 

1 . We use units GeV, TeV, PeV and EeV, increasing 
energy in steps of one thousand. 
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Figure 1 : The neutrino sky from the lowest energy neutrinos produced in the big bang to the highest energies 
associated with the sources of the cosmic rays, here assumed to be gamma ray bursts or, ahernatively, active 
galaxies. These will be the target of kilometer-scale neutrino detectors such as IceCube and KM3NeT. 
Neutrinos at intermediate energies, produced in the sun, supernovae and in collisions of cosmic rays in the 
atmosphere, have been studied by SuperK and similar detectors[8] 



neutrinos in highly transparent and well shielded 
deep water or ice. At the higher energies the neu- 
trino cross section grows and secondary muons 
travel up to tens of kilometers to reach the detector 
from interactions outside the instrumented volume; 
see Fig. 2. 

The construction of kilometer-scale instruments 
such as IceCube at the South Pole and the fu- 
ture KM3NeT detector in the Mediterranean, have 
been made possible by development efforts that 
resulted in the commissioning of prototypes that 
are two orders of magnitude smaller, AMANDA 
and ANTARES[5 |. Their successful technologies 



have, in turn, relied on pioneering efforts by the 
DUMAND[6] and Baikal[7], as well as the Macro 
and SuperK collaborations[8]. While much larger 
than the latter, kilometer scale neutrino telescopes 
are insensitive to neutrinos in the MeV-GeV en- 
ergy range and have a typical threshold of tens 
of GeV; this is the price one pays for reaching 
large detection volume. IceCube[9] is under con- 
struction and taking data with a partial array of 
1320 ten inch photomultipliers positioned between 
1500 and 2500 meter and deployed as beads on 22 
strings below the geographic South Pole. Its ef- 
fective telescope area already exceeds that of its 
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Figure 2: Neutrino telescopes take advantage of 
the large cross section of high-energy neutrinos 
and the long range of the muons produced. The 
detector consists of a lattice of photomultipliers de- 
ployed in a shielded and optically clear medium 
that is transformed into a Cherenkov detector. 

predecessor AMANDA by roughly one order of 
magnitude. The detector will grow by another 
14~18 strings in the 2007-08 Antarctic summer to 
be completed in 2011 with 80 strings [10]. 

The AMANDA experiment has observed neutrinos 
with energies as high as ^lOOTeV, at a rate con- 
sistent with the flux of atmospheric neutrinos ex- 
trapolated from lower energy measurements; see 
Fig. 1. The fluxes of cosmic neutrinos shown 
in the figure at higher energies are, in contrast, 
a matter of speculation. It is known that non- 
thermal sources such as supernova remnants, AGN 
and GRB accelerate electrons to energies close to 
1 Pe V. Their existence is inferred from the observa- 
tion of TeV gamma rays whose spectrum extends 
to ^ 100 TeV in some sources[ll]. All observa- 
tions can be accommodated in models where the 
origin of the photons, from radio to TeV energy, 
is synchrotron radiation by the electrons and, at 
the highest energies, inverse Compton scattering of 
ambient light, primarily the synchrotron photons 
themselves. There is no reason why non-thermal 
sources would not accelerate protons or nuclei 
along with the electrons, turning them into sources 
of cosmic rays; unfortunately, at this time, no evi- 
dence for such cosmic ray accelerators exists. A 
conclusive signature for the presence of cosmic 
rays in the sources is the production of pions on 



ambient radiation and matter. Pion production is 
revealed by observing the decay products, photons 
and neutrinos. While it has been a challenge to dis- 
entangle such pionic photons from those produced 
by purely electromagnetic processes [12], charged 
pions decaying into neutrinos yield incontrovert- 
ible evidence. The anticipated neutrino fluxes are 
shown in Fig. 1, assuming that AGN or, alterna- 
tively, GRB, happen to be the correct guess for the 
unknown sources of the cosmic rays. If not, the 
real sources may be revealed by neutrinos that, un- 
like charged primaries, point back to their site of 
origin. Neutrino astronomy must succeed because, 
after all, cosmic rays exist. The critical question 
is whether our estimate of the level of the neutrino 
fluxes associated with the observed cosmic rays is 
robust; it sets the scale of the detector. 

Thus the faith of neutrino astronomy is intertwined 
with cosmic ray physics beyond the traditional sub- 
ject of GZK neutrinos to which we will return later. 
While kilometer-scale neutrino detectors are dis- 
covery experiments with missions as diverse as 
particle physics and the search for dark matter - 
see Table 1 - their size as astronomical telescopes 
is very much anchored to the observed fluxes of 
galactic and extragalactic cosmic rays. Cosmic ac- 
celerators produce particles with energies in excess 
of 10^ TeV; we still do not know where or how. 
The flux of cosmic rays observed at Earth follows 
a broken power law; see Fig. 3. The two power 
laws are separated by a feature dubbed the "knee". 
Circumstantial evidence exists that cosmic rays, up 
to perhaps EeV energy, originate in galactic super- 
nova remnants. Any association with our Galaxy 
disappears in the vicinity of a second feature in the 
spectrum referred to as the "ankle". Above the an- 
kle, the gyroradius of a proton in the galactic mag- 
netic field exceeds the size of the Galaxy and it is 
generally assumed that we are witnessing the on- 
set of an extragalactic component in the spectrum 
that extends to energies beyond 100 EeV. Obser- 
vation of the GZK feature in the HiRes and Auger 
spectra near the energy threshold for pion produc- 
tion on microwave photons, provides further sup- 
port for the existence of an extragalactic compo- 
nent.^ While the enigmatic nature of the highest 

2. That the cutoff is associated with the upper energy 
range of the accelerator(s) can at this point not be ruled 
out. 
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energy cosmic rays has been widely advertised, it 
is also a fact that the origin of the galactic cosmic 
rays has not been established. 

Neutrinos Associated with Extragalac- 
tic Cosmic Rays 

It is routinely emphasized how small cosmic ray 
particle fluxes are; this may be besides the point. 
The energetics of the accelerator is likely to be 
more revealing. A hint in this direction comes 
from conventional astronomy. While the diffuse 
universal flux of photons falls by eighteen orders 
of magnitude between microwave and GeV-energy, 
the energy carried by the particle flux drops by less 
than five. Sources are known that emit most of 
their energy in TeV photons. The energy is the key. 
The argument has been well advertised for galac- 
tic cosmic rays where energetics points at their su- 
pernova origin. By integrating the observed flux 
in Fig. 3 we can obtain the energy density ps of 
cosmic rays in the galaxy from the relation that 
flux = velocity x density, or 
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The answer is that ^ 10~^^ ergcm"^. This 
is also the value of the corresponding energy den- 
sity B'^/Stt of the microgauss magnetic field in the 
galaxy. The accelerating power needed to maintain 
this energy density is 10~^^ erg/cm^^s given that the 
average containment time of the cosmic rays in our 
galaxy is 3 x 10^ years. For a nominal volume of 
the galactic disk of 10^^ cm^ this requires an ac- 
celerator delivering lO**^ erg/s. This happens to be 
10% of the power produced by supernovae releas- 
ing 10^^ erg, or the one percent of a solar mass that 
is not released in neutrinos, every 30 years. This 
coincidence is the basis for the idea that shocks 
produced by supernovae expanding into the inter- 
steUar medium are the origin of the galactic cosmic 
rays[14]. 

Let's follow the same logic for the extragalactic 
component in Fig. 3. The flux above the ankle 
is often summarized as "one 10^^ eV particle per 
kilometer square per year per steradian". This can 



be translated into an energy flux 



dE J (1010 cm2)(3 x 10^ sec) sr 

= 3 X lO-^GeVcm-^s-^sr-^ 

(2) 

FoUowing the procedure applied above to the 
galactic component we obtain an energy density of 
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taking the extreme energies of the accelerator(s) to 

be i^max/i^min ^ 10^. 

The energy content derived "professionally" by in- 
tegrating the spectrum in Fig. 3 assuming an E~'^ 
energy spectrum, typical of shock acceleration, 
with a GZK cutoff, is ~ 3 x 10" erg cm"^. This 
is within a factor of our back-of-the-envelope es- 
timate recalling that 1 TeV = 1.6 erg. The energy 
density represents the universe's filling factor in 
cosmic rays, equivalent to 410 microwave-energy 
photons per cubic centimeter. 

The power required for a population of sources to 
generate this energy density over the Hubble time 
of 10 1° years is ~ 3 x 10^ erg s~i per (Mpc)^ or, 
as often quoted in the literature, ~ 5 x 10*^ TeV 
per (Mpc)^ per year. This works out to[15] 

• 3 X 10^® erg s~i per galaxy, 

• 3 X 10^^ erg per cluster of galaxies, 

• 2 X 10^^ erg s~i per active galaxy, or 

• 2 X 10^1 erg per cosmological ganmia ray 
burst. 

The coincidence between the last two numbers 
with the observed output in electromagnetic ra- 
diation of these sources, explains why AGN and 
GRB emerged as leading candidates for the cos- 
mic accelerators. In either case, it suffices that the 
shocks associated with acceleration near the black 
hole dump roughly equal energy in electrons and 
protons to accommodate the observed coincidence, 
with the electron energy observed as radiation by 
synchrotron and inverse Compton scattering. 

For GRB the argument is reminiscent of the 
one favoring the "evidence" for galactic super- 
nova as cosmic ray accelerators. Observations 
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Table 1: Built as discovery instruments, neutrinos telescopes nevertheless target a range of particle and 



astrophysics problems. 




Favorite Sources 


Possible Science 


Atmospheric 

(~ 100, 000 per year, up to 1000 TeV, charm?) 


Oscillations 
New interactions 

Test of relativity and equivalence principle 


GRB 

(successful and failed) 


Sources of cosmic rays 

Test of Lorentz invariance 

Planck scale physics, quantum decoherence 


AGN 


Sources of cosmic rays 


Starburst galaxies 




Siiinpmova rpmnants 

Also microquasars, magetars, PWNe, binaries 
unidentified Egret sources, plane of the galaxy 


Soiirops of oosmio ravs 


Cosmic rays ineracting with microwave photons 


Identify sources of cosmic rays 
Neutrino cross section at EeV energy 


Dark matter 


Aimihilation in the sun, mostly spin-dependent 


Cosmic rays interacting with the sun 


Backrounds to WIMP search 


Supemovae explosion 


Deleptonization, TeV emission, hierarchy, sin ^13 



show, within one gigaparsec cubed, 300 GRB 
dumping about 10^^ erg of energy into the uni- 
verse in a single year. They therefore supply 
roughly 10*^erg/yr/Mpc^ in radiation and, as- 
suming equal energy in protons, we conclude that 
they represent an environment that can accommo- 
date the observed energetics of the extragalactic 
cosmic rays. A problem is that the same argument 
can be made to validate AGN as the sources of the 
highest energy cosmic rays; see above. In the end, 
the answer may lay elsewhere. At this point we 
should emphasize, in either case, the challenge re- 
mains to explain the acceleration of particles with 
energies as high as 10^ TeV, an energy which is in 
either source near the dimensionally allowed upper 
Umit set by the Hillas formula. 

Where do neutrinos fit into this? The assump- 
tion that the energy in neutrinos coincides with the 
matching energies observed in electromagnetic ra- 
diation and cosmic rays yields the level of neutrino 
fluxes associated with the cosmic rays shown in 
Fig. 1. In the end the neutrino flux is therefore 



the flux of Eq. 2.^ It is often referred to as the 
Waxman-BahcaU "bound"[16]. A source creating 
equal fluxes of cosmic rays, gamma rays and neu- 
trinos is rather generic; see Fig. 4. The usual as- 
sumption is that cosmic rays are accelerated in a 
region of high magnetic fields, most likely associ- 
ated with shocked particle flows driven by the grav- 
ity of a black hole. They interact with the radiation 
fields surrounding the black hole. The most impor- 
tant processes are p -|- 7 A+ ^ tt'^ + p and 
p -h 7 ^ A+ TT"*" + n. While the secondary 
protons may remain trapped in the acceleration re- 
gion, roughly equal numbers of neutrons and de- 
cay products of neutral and charged pions escape. 
The energy escaping the source is therefore dis- 
tributed among cosmic rays, gamma rays and neu- 

3. There are many corrections to the equality, from the 
details of the particle physics, including neutrino oscil- 
lations, to the fact that cosmic rays only reach us from 
within a GZK absorption length of ~ 50 Mpc while neu- 
trinos travel unimpeded from sources at all redshifts. In 
the end, these modifications cancel "within a factor". 
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Figure 3: At the energies of interest here, the cosmic ray spectrum consists of a sequence of 3 power laws. 
The first two are separated by the "knee" (left panel), the second and third by the "ankle". There is evidence 
that the cosmic rays beyond the ankle are a new population of particles produced in extragalactic sources; 
see right panel. 



trinos produced by the decay of neutrons, neutral 
pions and charged pions, respectively. This generic 
scenario accommodates the observation of equal 
energy in cosmic rays and electromagnetic radia- 
tion, and extends it to neutrinos. Clearly both GRB 
and AGN environments can accommodate this sce- 
nario although with very dissimilar black holes and 
radiation targets for pion production. If we take 
this picture seriously, our previous estimate must 
be corrected for the fact that the pion takes only 
25% of the energy of the secondary neutron thus 
changing the energy balance between cosmic rays 
and neutrinos and reducing their flux. In the end 
we estimate that the muon-neutrino flux associated 
with the sources of the highest energy cosmic rays 
is loosely confined to the range 

(3) 

depending on the cosmological evolution of the 
cosmic ray sources. Model calculations assuming 
that active galaxies or gamma-ray bursts are the ac- 
tual sources of cosmic rays yield event rates simi- 



lar to the generic energetics estimate presented; see 
Fig. 5. 

The anticipated neutrino flux thus obtained 
has to be compared with the limit of 7.4 x 
10^® GeV cm^^s^^sr^^ reached after the first 
4 years of operation of the completed AMANDA 
detector in 2000-2003 [17]. On the other hand, 
for conservative assumptions for the charm back- 
ground and for the detector performance, IceCube 
has the capability to observe a flux that is an order 
of magnitude below this limit with 5 cr significance 
in 3 years [9]. The exact value of the IceCube sen- 
sitivity depends on the magnitude of the dominant 
high energy neutrino background from the prompt 
decay of atmospheric charmed particles. The level 
of this background is difficult to anticipate theoret- 
ically and little accelerator data is available in the 
energy and Feynman-x range of interest[18]. 

The neutrino event rate is obtained by folding the 
flux predicted with the probability that the neutrino 
is actually detected in a high energy neutrino tele- 
scope; only one in a million neutrinos of TeV en- 
ergy interacts and produces a muon that reaches 
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Figure 5: Our estimate of the flux of neutrinos associated with the sources of the highest energy cosmic rays 
(the shaded range labeled WB) is compared to the limits established by the AMANDA experiment reached 
with 800 days of data[17]. AMANDA'S sensitivity is within a factor of 2 of the most optimistic predictions. 
Also shown are fluxes predicted by specific models of cosmic ray accelerators: active galaxies labeled 
StSa[19] and MPR[20], GRB[21] and the diffuse flux produced by cosmic ray producing active galaxies 
on microwave photons[22] labelled RB. Data for the background atmospheric neutrino flux are from the 
AMANDA experiment. The IceCube experiment will be sensitive to all predictions after a few years of 
operation of the full detector. It has sensitivity to the larger fluxes by operating the partially completed 
detector that already now exceeds AMANDA in instrumented volume. 



the detector. This probability is given by the ratio 
of the muon and neutrino interaction lengths in the 
detector medium, Ap/Ai/[4] and therefore grows 
with energy. For the flux range estimated above we 
anticipate 100-500 detected muon neutrinos per 
km^ per year, with the higher range close to what 
is already ruled out by AMANDA. In any case, the 
lower value represents the more realistic estimate 
as previously argued and the 100 events predicted 
will be further reduced by the realities of reject- 
ing detector backgrounds, especially for spectra 
steeper than the E^'^ assumed throughout. On the 
other hand, IceCube's effective area for muon neu- 



trinos exceeds 1 km^ and equal fluxes of electron 
and tau neutrinos are expected[9]. If IceCube con- 
struction remains on schedule, the instrument will 
accumulate 1 km^ year of data within the next two 
years; the confrontation of these arguments with 
data is imminent, certainly on the 40 year timescale 
it took to develop the technology for the detectors. 

We next return to galactic cosmic rays. Since 
the last cosmic ray conference, new observations 
by air Cherenkov telescopes as well as the results 
from an all-sky survey by the Milagro detector, 
are suggestively pinpointing supernova remnants 
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Figure 4: Cosmic beam dumps exist: sketch of 
cosmic ray accelerator producing photons. The 
charged pions that are inevitably produced along 
with the neutral pions will decay into neutrinos. 

as the sources. Some would say that the smoking 
gun is still missing and neutrinos may be the key. 

Cosmic Neutrinos Associated with 
Galactic Supernova Remnants 

Can kilometer-scale neutrino detectors observe 
neutrinos pointing back at the accelerators of the 
galactic cosmic rays? It is believed that galac- 
tic accelerators are powered by the conversion of 
10^" erg of energy into particle acceleration by dif- 
fusive shocks associated with young (1000-10,000 
year old) supernova remnants expanding into the 
interstellar medium. The cosmic rays will interact 
with hydrogen atoms in the interstellar medium to 
produce pions that decay into photons and neutri- 
nos. These provide us with indirect evidence for 
cosmic ray acceleration. The new twist here is 
that the eventual observation of pionic gamma rays 
allows for a straightforward determination of the 
neutrino flux. 

The HESS telescope has opened a new era in as- 
tronomy by producing the first resolved images of 



sources in TeV gamma rays, particularly, in this 
context, of the supernova remnant RX J1713.7- 
3946 [23]. While the resolved image of the source 
reveals TeV gamma ray emission from the whole 
supernova remnant, it shows a clear increase of the 
flux in the directions of known molecular clouds. 
This is suggestive of protons, shock-accelerated in 
the supernova remnant, interacting with the dense 
clouds to produce neutral pions that are the ori- 
gin of the observed increase of the TeV photon 
signal. The magnitude of the photon flux is con- 
sistent with a site where protons are accelerated 
to energies typical of the main component of the 
galactic cosmic rays. A similar extended source 
of TeV gamma rays tracing the density of molecu- 
lar clouds has been identified near the galactic cen- 
ter Protons, apparently accelerated by the remnant 
HESS J1745-290, diffuse through nearby molecu- 
lar clouds to produce a signal of TeV gamma rays 
that traces the density of the clouds [24]. Fitting 
the observed spectrum by purely electromagnetic 
processes is challenging because the relative height 
of the inverse Compton and synchrotron peaks re- 
quires very low values of the S-field, inconsis- 
tent with those required to accelerate the electron 
beam to energies that can accommodate the obser- 
vation of 100 TeV photons. Nevertheless, an ex- 
clusively electromagnetic explanation of the non- 
thermal spectrum is not impossible, even favored 
by some [25]. One can, for instance, partition the 
remnant in regions of high and low magnetic fields 
that are the respective sites of acceleration and in- 
verse Compton scattering. 

Supernovae associated with molecular clouds are 
a common feature of associations of OB stars that 
exist throughout the galactic plane. Although not 
visible to HESS, possible evidence has been accu- 
mulating for the production of cosmic rays in the 
Cygnus region of the galactic plane from a vari- 
ety of experiments [28, 31, 32, 33, 34]. Most in- 
triguing is a Milagro report of an excess of events 
from the Cygnus region at the 10.9 a level [34]. 
The observed flux within a 3° x 3° window is 
70% of the Crab at the median detected energy 
of 12 TeV and is centered on a source previ- 
ously sighted by HEGRA. Such a flux largely ex- 
ceeds the one reported by the HEGRA collab- 
oration, implying that there could be a popula- 
tion of unresolved TeV 7-ray sources within the 
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Cygnus OB2 association. In fact, they report a 
hot spot, christened MGRO J2019+37, at right as- 
cension = 304.83° ± 0.14stat ± 0.3sys and dec- 
Unation = 36.83° ± 0.08stat ± 0.25sys [34]. A 
fit to a circular 2-dimensional Gaussian yields a 
width of 0.32 ± 0.12 degrees, which for a dis- 
tance of 1.7 kpc suggests a source radius of about 
9 pc. As the brightest hotspot in the Milagro map 
of the Cygnus region, it represents a flux of 0.5 
Crab above 12.5 TeV. A model proposed [35] for 
MGRO J2019H-37 is that of a cosmic ray beam 
which escapes from the OB star cluster and inter- 
acts with a molecular cloud positioned a few de- 
grees to the southeast. Interestingly, the Tibet AS- 
gamma Collaboration has observed a cosmic ray 
anisotropy from the direction of Cygnus, which is 
consistent with Milagro's measurements [36]. As 
for the HESS sources, these observations suggest 
the production of cosmic rays as well as a variety 
of opportunities for neutrino production. 

If the TeV gamma ray signals are indeed of pi- 

onic origin, only particle physics establishes the 
rate of the accompanying neutrinos. Proton-proton 
collisions yield two charged pions for every neu- 
tral pion, with every charged pion decaying into a 
muon neutrino and antineutrino (one from the pion 
decay, the other from the decay of the secondary 
muon) and a neutral pion into two photons. So, the 
muon neutrino flux would be equal to the photon 
flux were it not for a factor two reduction from os- 
cillations. The prediction is simple; to first order 
there is one muon neutrino for every photon pro- 
duced in the source. Because the protons transfer 
on average 20% of their energy to secondary pi- 
ons, and the four leptons in the charged pion decay 
chain tt — > > e + + f^) + take roughly 
equal energy, neutrinos with 0.05 of the cosmic ray 
energy are produced. Similarly, photons with 10% 
of the proton energy originate from the decay of 
neutral pions. Accelerators producing cosmic rays 
reaching the "knee" must produce photons with en- 
ergies up to 100 TeV and neutrinos up to half that 
energy. This requirement is consistent with obser- 
vations of RX J 17 1 3 .7-3946 and MGRO J20 1 9h-37 
discussed above. They are the targets for neutrino 
observation of neutrino telescopes located in the 
Southern and Northem hemispheres, respectively. 

Whereas the relation between neutrino and gamma 
ray fluxes is direct, the information on their spec- 



trum is often limited. This is especiaUy the case 
for the hotspot MGRO J2019H-37 where we have 
to model the spectrum on the basis of a measure- 
ment at a single energy; the spectral slope has not 
been measured. Uncertainties in the calculation are 
associated with the propagation of the cosmic rays, 
with the value of the magnetic fields, and the age 
of the remnant. After investigating the wide pa- 
rameter space of models for MGRO J2019H-37, it 
has been shown[35] that the neutrino flux can nev- 
ertheless be predicted within a factor of 2 once the 
model flux is normalized at 12.5 TeV to the Mila- 
gro data and a Umit at GeV energy is imposed re- 
flecting the fact fliat EGRET did not observe a GeV 
counterpart [37]. The range of neutrino fluxes and 
event rates in IceCube are shown in Figure 7 as- 
suming a 2.2 injected slope. The rates are within 
the range 2 < dN/dt < 3.8 events per year 
with the IceCube threshold at 50 GeV. The range 
is bound by the fact that the Milagro observation 
strongly constrains the flux in the energy range 1- 
20 TeV, where the neutrino detection probability is 
highest, resulting in similar predictions for dissim- 
ilar SNR characteristics. 

The irreducible atmospheric background, due to 
neutrinos produced in the Northern atmosphere in 
cosmic ray showers, is calculated using the results 
of Ref. [38]. In 15 years of operation we predict 
4.5 a < N/^/N^tmo < 7.7 a and if the higher 
end of the predicted event rate range is realized, 
5 (T is possible in 4.3 years. 

We note that the Milagro collaboration [39] has re- 
cently detected multiple additional sources besides 
MGRO J2019H-37, most with fluxes close to 0.5 
Crab. The sources with possible counterparts in the 
GeV range indicate a spectral index of ~ —2.3. If 
we compute the flux of neutrinos from the Milagro 
sources (not including the Crab Nebula) detected 
with post-trial significance of greater than 5 a as- 
suming a power-law index of —2.3, we get a total 
event rate in IceCube of 6.9 neutrinos/year. If we 
also include the more tentative sources, the event 
rate increases to 11.5 neutrinos/year In the long 
run, a correlation analysis of the IceCube and Mi- 
lagro skymaps should make the detection of these 
sources Ukely. 

It is important to emphasize that the photon flux 
from the Milagro sources is consistent with the 
flux expected from a typical cosmic ray generat- 
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Figure 6: 7-ray spectra with injection s = —2.2. The black dashed Hne is for a magnetic field of 0.1 fiG 
and age of 1400 years, the red dotted line is for a magnetic field of 1 fiG and age of 5000 years, the green 
dash-dotted line is for a magnetic field of 10 /iG and age of 5000 years, and the blue solid line is for a 
magnetic field of 50 /iG and age of 5000 years. 



ing supernova remnant interacting with the inter- 
stellar medium[14]. In other words, the TeV flux 
is consistent with the energetics that are required 
to power the cosmic ray flux in the galaxy. Al- 
ternative candidates have been suggested for the 
sources of the galactic cosmic rays, for instance 
microquasars. Reversing the argument for super- 
nova remnants, cosmic ray energetics requires that 
they should have left their imprint on the Milagro 
skymap and they did not. It is very suggestive that 
the Milagro sources are the cosmic ray accelera- 
tors. 

Greisen-Zatsepin-Kuzmin Neutrinos 

From our previous discussions the case for a 
kilometer-scale neutrino detector clearly emerges 
even though the predicted fluxes are hardly guaran- 
teed. Neutrinos are guaranteed when both the ac- 
celerator beam and the pion producing target mate- 
rial can be identified. This is the case for GZK neu- 



trinos produced by extragalactic cosmic rays inter- 
acting with microwave photons. The event rate is 
of order one per year in a kilometer-cubed detector; 
even kilometer-scale detectors such as IceCube are 
marginal in this case and unlikely to accumulate a 
statistically useful sample of events. It is also im- 
portant to be aware of the fact that this event rate 
is only determined within large ambiguities asso- 
ciated with the calculation of the flux and with the 
efficiency of the detectors at such high energies, 
typically EeV. 

A suite of experiments have been searching for 
these rare events exploiting the fact that the 
showers initiated by EeV neutrinos emit coherent 
Cherenkov radiation in the 20AIHz ^ IGHz ra- 
dio range[40]. EeV neutrinos can also be detected 
as near-horizon airshowers. The large pathlength 
in the atmosphere of near-horizon cosmic rays 
gives Auger the capabilities of an underground 
experiment; interesting limits have already been 
obtained[41]. The possibility of exploiting this en- 
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s=-2.2, T=1400 years, B=0.1(iG 




Figure 7: Events due to the gamma-ray spectra with injection a = —2.2 shown in Fig. 6. 



ergy range with acoustic detection techniques is 
also under intense investigation[43]. 
As early as 1962 Askaryan proposed using ice 
as a Cherenkov medium for neutrino-induced ra- 
diation of GHz wavelengths [42]. In the early 
1990s we witnessed a renewed interest in Askar- 
ian's proposal with the recognition that the rela- 
tively high neutrino energy threshold of the tech- 
nique, 10 PeV or more in a reasonably scaled em- 
bedded detector in ice and even higher for other ge- 
ometries, is well-matched to a number of physics 
goals, most notably the detection of GZK neutri- 
nos. The RICE collaboration pioneered the tech- 
nique by positioning dipole antennas in AMANDA 
holes at depths of several 100 m [44]. A suite 
of experiments followed, GLUE and FORTE [45], 
setting the first limits at extremely high ener- 
gies exceeding 10^" eV. More recent detectors in- 
clude Lunaska[46] and the ANITA balloon pay- 



load, which completed a prototype flight in 2004, 
and its first full-payload flight in early 2007. 

Because of absorption by the earth, neutrinos of 
EeV energy predominantly produce signals in the 
ice near the horizon. These can refract and be de- 
tected at the ice surface or, as is the case for the 
ANITA experiment, by an array of 32 quad-ridged 
horn antennas floating in the circumpolar wind 
from a balloon 30 km above Antarctica. It scans 
a volume of order 10^ km"^ of radio-transparent ice 
for neutrino interactions over a continent that pro- 
duces very little radio noise, unlike populated re- 
gions of the world. Data from the first flight with 
the completed detector are subject to a blind analy- 
sis from which the results are eagerly awaited[47]. 

The emission of GHz Cherenkov radiation from 
electromagnetic cascades is counterintuitive be- 
cause, according to the Frank-Tamm formula, the 
power in the signal should be suppressed by a fac- 
tor proportional to the frequency of 10^ relative to 
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light, and is therefore unobservable. Askarian re- 
alized that this argument is wrong. At MHz~GHz 
frequencies the wavelength sampled exceeds the 
dimensions of the shower and the radiation is 
therefore coherent. The power coherently radi- 
ated is proportional to the square of the number of 
shower particles times their charge, Ne. Coherence 
therefore compensates for the suppression in fre- 
quency provided the number of shower particles is 
sufficiently large; we now know that it is detectable 
above the thermal background in ice provided the 
shower energy exceeds ~ 10 PeV[48]. There is an- 
other subtlety here because a shower, consisting of 
electron-positron pairs and photons, is to a first ap- 
proximation electrically neutral, Ne = after sum- 
ming over all particles. Askarian reahzed that an 
electromagnetic cascade develops a ~ 20% excess 
of negative charge, predominantly produced by the 
Compton scattering of the large number of MeV 
shower photons on atomic electrons. Although the 
complete calculation of the effect is significantly 
more complicated, his estimate of the charge ex- 
cess was quahtatively correct[48]. The existence 
of the coherent signal was demonstrated by ex- 
periments in the SLAC beam using sand and ice 
targets [49]. The radiated power and its energy de- 
pendence measured agreed with the modem calcu- 
lations. 

There are additional considerations that make this 
technique attractive. The absorption length of ra- 
dio waves in the cold ice is close to 1 km up 
to 1 GHz frequency [50J, to be contrasted with ~ 
100 m for light, and their detection with antennas 
relatively simple. Moreover, it may be possible to 
perform the experiment with antenna arrays posi- 
tioned at or near the surface avoiding the time con- 
suming and expensive necessity of drilling. The 
array would consist of stations that can be visu- 
alized as descoped versions of the ANITA exper- 
iment, simpler yet having the possibility to make 
a stand-alone detection of the shower direction 
and energy. Already in the early 1970s Gusev 
and Zheleznykh[51] proposed a surface radio array 
with a 10 km^ footprint to detect of order 10 PeV 
neutrinos via antennas with grid spacing of several 
hundred meters. 

The new Auger data presented at this meeting 
have cast GZK searches in a somewhat different 
Ught. The confirmation of a cutoff in the spec- 



trum imphes that the rate of super-lO^" eV cosmic 
ray events is low, of order one per year. Proton 
astronomy, while certainly possible, will be sta- 
tistically challenging. An experiment observing 
several GZK neutrinos per year is complementary 
and the neutrinos are guaranteed to point at their 
sources. One can indeed extract directional infor- 
mation from a neutrino that is produced within a 
GZK radius of ~ 50 Mpc of a source that is at a 
distance of hundreds to thousands of megaparcecs. 
On the negative side, measurements of the depth 
of shower maximum (ccmax) suggest the possibil- 
ity that a fraction of the primaries are heavy which 
will reduce the GZK neutrino flux. Heavy pri- 
maries photodisintegrate on the universal photon 
background before photoproducing the pions that 
are the parents of GZK neutrinos. 

In any case, we have to be realistic about the confi- 
dence with which we can anticipate the GZK flux. 
Despite the availability of significantly improved 
data on the spectrum near and beyond the "an- 
kle", these can still be accommodated with a wide 
range of assumptions regarding the injection slope 
at the source, the cosmological evolution of the 
sources and, as already mentioned, the composi- 
tion at injection. At this point scenarios yielding 
GZK fluxes ranging from imobservable to several 
per km^ year can be entertained[52], see Table 2 
taken from reference[52]. 

Independent of the difficulty of anticipating the 
flux, the GZK event rate is proportional to the neu- 
trino cross section whose calculation is challenging 
at these high energies where accelerator data pro- 
vide little guidance. That it is poorly known has 
been routinely ignored in the radio business and 
neutrino Umits from existing experiments can cer- 
tainly be weakened by arguing for a less optimistic 
extrapolation of the neutrino cross section. In the 
Standard Model the neutrino cross section is sim- 
ply proportional to the g(x, Q^) + q{x, Q"^) quark 
structure function, with Q"^ ~ ^w- quarks 
produced by gluons dominate the distribution func- 
tion because, at these high energies, the relevant 
values of the fractional momenta of the quarks be- 
come as small as x ~ 10~^. For the relevant 
range, HERA data barely constrain the par- 
ton distributions to x ~ 10~^, while collider data 
on inclusive jet production provide indirect evi- 
dence to 10~^. Perturbative QCD dictates a power 
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Table 2: The rates of muon and shower events initiated by GZK neutrinos in a kilometer- scale neutrino 
telescope (such as IceCube or KM3) for a range of choices of injected spectra and chemical composition 
consistent with both the Auger spectrum and Xmax measurements. For comparison, we also show the event 
rates for the case of an all-proton spectrum with injected spectral index 2.2 and £^max = lO^^ eV. Relative 
to the all-proton case, models are consistent with the Auger data which range from almost no difference, to 
approximately two order of magnitude suppression. 



Composition 


a 


E^^JZ (eV) 


Muons (km 2yr 1) 


Showers (km'^ yr— 1) 


100% N 


1.6-1.9 


10^2 


0.20-0.0081 


0.15-0.0064 


100% Si 


1.6-2.0 


1022 


0.21-0.045 


0.16-0.035 


100% Fe 


1.6-2.1 


1022 


0.11-0.014 


0.085-0.012 


100% Fe 


1.4-1.7 


1021 


0.019-0.0076 


0.017-0.0075 


50% N, 50% p 


1.8-2.1 


1022 


0.23-0.013 


0.18-0.10 


50% Si, 50% p 


1.6-2.1 


1022 


0.30-0.095 


0.220-0.075 


50% Si, 50% p 


1.4-1.5 


1021 


0.059-0.051 


0.050-0.043 


7% Si, 93% p 


2.0-2.2 


1022 


0.69-0.66 


0.52-0.50 


2% Si, 98% p 


1.4-1.8 


1021 


0.75-0.59 


0.60-0.47 


50% Fe, 50% p 


1.6-2.1 


1022 


0.15-0.043 


0.11-0.034 


10% Fe, 90% p 


1.4-1.9 


1021 


0.14-0.10 


0.11-0.080 


3% Fe, 97% p 


2.1 


1022 


0.68 


0.51 


1% Fe, 99% p 


1.4-1.9 


1021 


0.74-0.53 


0.59-0.43 


100% p (for comparison) 


2.2 


1022 


0.76 


0.60 



law behavior in x that allows us to extrapolate to 
smaller values. But what is its range of appUca- 
bility? Clearly the growth of the number of con- 
stituents built into the perturbative extrapolation, 
and routinely assumed, cannot continue because it 
eventually violates the Froissart bound on the cross 
section. Screening of gluons will prevent this from 
happening. There is no consensus on the energy 
where this screening of the gluon constituents in- 
side the proton, similar to that of nucleons inside 
nuclei, sets in, and what the mechanism is, possi- 
bly new physics such as the color glass condensate 
observed at RHIC. A logarithmic reduction of the 
cross section has been argued for on the basis of a 
reanalysis of the HERA data[53]. 

Given the challenges in determining the flux as 
well as the acceptance of the detectors, one should 
reasonably conclude that the rate of GZK events 
must be determined experimentally. Ideally, the 
initial data of the IceCube or ANITA experiments 
would give us an indication. Continued data tak- 
ing of IceCube and a planned ANITA flight in two 
years will establish GZK rates that cover a range 
of optimistic expectations in Table 2. For instance, 
each experiment would detect a few events cor- 



responding to the "reference" flux of Engel et al. 
which assumes a (1 + z)^ redshift evolution of the 
sources[54]. In the absence of information from 
the existing experiments a straightforward way to 
proceed is to upgrade the IceCube experiment us- 
ing the optical Cherenkov technique and the in- 
frastructure already available. The possibility has 
been investigated in reference [55] and the most 
straightforward conclusion is that significant up- 
grades are required, in fact on the scale of IceCube 
itself, to increase the event rate by "a factor". The 
necessity to switch from photomultipUers to radio 
antennas strongly suggests itself. 

The proposal is for staged deployments of radio de- 
tectors on an approximately kilometer grid expand- 
ing on the hexagonal outlay of the IceCube strings 
in stages of 6, 12,... radio detectors. If the event 
rate established with the initially deployed detec- 
tors warrants it, this array can be expanded into an 
instrument exceeding 100 km^ in effective volume 
detecting hundreds of events for doing cosmic ray 
and particle physics. It would open up the possibil- 
ity to measure the neutrino cross section at EeV en- 
ergy. This is of obvious interest to particle physics; 
the case has been extensively iUustrated in the con- 
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text of TeV- scale gravity. The approach has the 
critical advantage that a fraction of the events will 
also be detected by IceCube allowing for a calibra- 
tion of this novel detection method. As many as 
20% percent of the events would be hybrid events 
seen by both detectors provided one surrounds Ice- 
Cube by an additional a single hexagon of opti- 
cal strings at a distance of 500 m as suggested in 
reference[55]. 

Whereas calibration of these experiments would be 
desirable, there should no longer be any doubt that 
the radio technique, suggested by Askarian[42] al- 
most half a century ago, is robust. The signal pro- 
duced by neutrinos interacting in the ice can be 
calculated — it is just QED — and the theoretical 
results have been validated with accelerator exper- 
iments. There are of course additional hurdles to 
overcome in establishing the sensitivity of experi- 
ments that, unlike IceCube, have a threshold that is 
too high to allow for calibration using atmospheric 
neutrinos. It is important to realize that the consid- 
erable problems to establish the radio technique for 
detecting air showers are not relevant here, for in- 
stance the effects of the earth magnetic field which 
are much more difficult to quantify, are negligible. 

The ARIANNA project[56] proposes a variation 
on the concept for a surface radio array by position- 
ing detector stations on the Ross Ice Shelf where 
the Antarctic continent is supported by water rather 
than rock. The experiment will capitalize on the 
fact that at radio frequencies the water-ice bound- 
ary below the shelf act as a mirror reflecting sig- 
nals produced by neutrino signals produced in any 
downward direction. The concept consists of sta- 
tions of cross-polarized antennas facing downward 
just below the snow surface positioned on a 100 by 
100 grid with 300 m spacing. 

Looking Forward 

While neutrino "telescopes" are discovery instru- 
ments with a variety of missions, the hope is that 
they may contribute to the resolution of the century 
old puzzle of the origin of cosmic rays, either by 
the detection of GZK neutrinos or by directly ob- 
serving neutrinos from the accelerators. Between 
now and the next International Cosmic Ray Con- 
ference we can look forward to 



• The results of the first ANITA flight with a 
complete pay load of antennas. 

• Results from IceCube that, by operating the 
detector as it grows, will reach a km?' year 
aperture. 

• Enhanced sensitivity of the Auger experi- 
ment to neutrinos that initiate horizontal air 
showers. 

• The initial design of a kilometer- scale neu- 
trino detector in the northern hemisphere. 

• Data as well as a wealth of ideas and initia- 
tives on detecting the acoustic and radio sig- 
natures produced by GZK neutrinos in ice, 
water, salt, permafrost . . . 
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